Prepared for submission to JCAP 



Global fits of the cMSSM including 
the first LHC and XENON100 data 



^ ; Gianfranco Bertone" David G. Cerdeno 5 Mattia Fornasa cl 

Roberto Ruiz de Austri d Charlotte Strege e Roberto Trotta 6 >f 

(N 



c3 



^ 



in 



"GRAPPA Institute, University of Amsterdam, Science Park 904, 1090 GL Amsterdam, 

Netherlands 

Departamento de Fisica Teorica and Instituto de Fisica Teorica UAM-CSIC, Universidad 

0^ \ Autonoma de Madrid, Cantoblanco, E-28049 Madrid, Spain 

— —t 

c Instituto de Astrofisica de Andalucia (CSIC), E-18008, Granada, Spain 

Instituto de Fisica Corpuscular, IFIC-UV/CSIC, Valencia, Spain 

Q_i' e Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Consort Road, 

£^. London SW7 2AZ, UK 

■'African Institute for Mathematical Sciences, 6 Melrose Rd, Muizenberg, 7945, Cape Town, 

i__i| South Africa 

cV 

Abstract. We present updated global fits of the constrained Minimal Supersymmetric 
Standard Model (cMSSM), including the most recent constraints from the ATLAS and CMS 
C" detectors at the LHC, as well as the most recent results of the XENON100 experiment. Our 

robust analysis takes into account both astrophysical and hadronic uncertainties that enter 
in the calculation of the rate of WIMP-induced recoils in direct detection experiment. We 
study the consequences for neutralino Dark Matter, and show that current direct detection 
data already allow to robustly rule out the so-called Focus Point region, therefore demon- 
strating the importance of particle astrophysics experiments in constraining extensions of 
the Standard Model of Particle Physics. We also observe an increased compatibility between 
results obtained from a Bayesian and a Frequentist statistical perspective. We find that 
upcoming ton-scale direct detection experiments will probe essentially the entire currently 
favoured region (at the 99% level), almost independently of the statistical approach used. 
Prospects for indirect detection of the cMSSM are further reduced. 
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1 Introduction 

The constrained Minimal Supersymmetric Standard Model (cMSSM) [1, 2] is the most widely 
discussed extension of the Standard Model (SM) of particle physics. Despite its relative 
simplicity, this model has the advantage of capturing some key phenomenological feature 
of Supersymmetry (SUSY), while making definite predictions for the properties of the the 
lightest neutralino Xi (which we denote x here, for simplicity) , a linear superposition of the 
superpartners of the neutral Gauge bosons and the neutral Higgses, which is by far the most 
popular Dark Matter (DM) candidate [3-6]. 

Recently, the ATLAS and CMS experimental collaborations have published the first 
search for supersymmetry (SUSY) [7-11], based on i/s = 7 TeV proton-proton collisions at 
the LHC recorded during 2010 (see the experimental papers above for details on the event 
topologies, selection cuts, etc.). Since no excess above the SM predictions has been observed, 
these searches have set new interesting constraints on physics beyond the SM, and several 
groups of authors have already studied the impact of these results on the cMSSM and more 
in general on SUSY (see e.g. Refs. [12-19]). 

Here, we perform new global fits of the cMSSM, following the methodology outlined 
in Refs. [20-22], where LEP constraints, precision tests of the SM and the most recent 
cosmological constraint on the relic abundance of the DM were implemented using nested 
sampling as a scanning algorithm, in order to derive both the most probable regions of 
the cMSSM (Bayesian) and the best fit regions (frequentist). We update the analysis to 
include new -B-physics observables, as well as both the aforementioned LHC constraints 
and the upper limit on the WIMP-nucleon scattering cross-section arising from DM direct 
detection experiments (see e.g. Ref. [23] and references therein). We implement the findings 



of the XENON100 collaboration [24], which recently reported the results of a two-phase time 
projection chamber with 48 kg of liquid xenon as ultra-low background fiducial target. In 
100.9 live days of data, acquired between January and June 2010, the collaboration recorded 
three candidate events, with an expected background of 1.8 =L 0.6 events. The rate of DM- 
induced recoil events in a direct detection experiment being directly proportional to the 
spin- independent scattering cross section of neutralinos off protons o"yL p , the upper limit on 
the recoil rate translates into an upper limit on a^_ „. 

The predictions for c5_ p within the cMSSM (given constraints other than direct detec- 
tion data) span approximately the range between 10 -7 pb to 10 -10 pb, and since XENON100 
excludes cr^_ p larger than 7.0 x 10~ 9 pb for a DM particle mass of 50 GeV at the 90% confi- 
dence level, one sees immediately that this has the potential to rule out significant portions 
of the cMSSM, under the hypothesis of a standard expansion rate of the Universe, and that 
the neutralinos make up most of the DM in the Universe. 

Particular care must be paid, however, when translating the information provided by a 
direct detection experiment to constraints on the cMSSM parameter space, since this requires 
specific modeling of the velocity distribution of DM particles and of the DM density in the 
solar neighborhood. As a consequence, the upper limit from XENON100 can be considered 
robust only if the uncertainties on these quantities are taken into account. Moreover, the value 
of OyLp depends on nuclear physics quantities, namely the hadronic matrix elements, which 
parametrize the quark composition of the proton, and this affects the way direct detection 
data are applied to constrain cMSSM parameters. Therefore in the present analysis we 
incorporate the local density, the quantities defining the DM velocity distribution and those 
modeling the proton composition as nuisance parameters, which are then either marginalized 
(i.e., integrated out) or maximised over to obtain inferences on the cMSSM parameters that 
fully account for those astrophysical and hadronic uncertainties. 

The paper is organized as follows: in Section 2 we present our theoretical framework, in 
Section 3 we present the formalism for the computation of recoil events, our modelization of 
the DM velocity distribution and our implementation of the XENON100 data. Section 4 is 
devoted to the presentation of the results, while prospects for the discovery of SUSY in the 
cMSSM are discussed in Section 5. We conclude in Section 6. 

2 Theoretical and statistical framework 

2.1 Model, priors and scanning algorithm 

The five parameters defining a generic cMSSM are i) the coefficient of the universal mass 
term for scalars mo and it) for gauginos mift, Hi) the coefficient of the trilinear interaction 
Aq, iv) the ratio between the vacuum expectation values of the Higgs bosons tan/3 and v) 
the sign of the mass term for Higges sgn(/x). We consider only four of these in our scan, 
fixing the sign of the \i coefficient to be positive, a choice favoured by the results on the 
anomalous magnetic moment of the muon. The scanned range of the remaining four cMSSM 
parameters is summarized in Table 1. Moreover, since residual experimental uncertainties 
on the value of some SM quantities can have strong effects on the values for the observables 
we are interested in [25], we also include in our scan the four SM parameters indicated in 
Table 2, which are considered as nuisance parameters. 

The cMSSM parameters together with the nuisance parameters in our scan are indicated 
by 0. The core of Bayesian statistics is Bayes' theorem, which reads (see Ref. [32] for further 



cMSSM Parameters 

Flat priors Log priors 

m [GeV] (50.0, 4000.0) (10 1 - 7 , 10 3 - e ) 
mi/2 [GeV] (50.0, 4000.0) (10 L7 , 10 3 - 6 ) 
A [GeV] (-4000.0, 4000.0) 

tan/3 (2.0, 65.0) 



Table 1. cMSSM parameters and the range covered on the scan. Flat priors are uniform in the masses, 
while log priors are uniform in the logarithm of the masses. 



SM nuisance parameters 




Gaussian prior 


Range scanned 


Ref. 


M t [GeV] 


173.1 ±1.3 


(167.0, 178.2) 


[26] 


m b {m b ) MS [GeV] 


4.20 ±0.07 


(3.92, 4.48) 


[27] 


[a em {M z ) MS ]- 1 


127.955 ± 0.030 


(127.835, 128.075) 


[27] 


a s (M z ) MS 


0.1176 ±0.0020 


(0.1096, 0.1256) 


[28] 


Astrophysical nuisance 


parameters 




P\oc [GeV/cm 3 ] 


0.4 ±0.1 


(0.001, 0.900) 


[29] 


fl sr [km/s] 


230.0 ± 30.0 


(80.0, 380.0) 


[29] 


v csc [km/s] 


544.0 ± 33.0 


(379.0, 709.0) 


[29] 


v d [km/s] 


282.0 ±37.0 


(98.0, 465.0) 


[29] 


Hadronic nuisance parameters 


fTu 


0.02698 ± 0.002 


(0.010, 0.045) 


[30] 


frd 


0.03906 ± 0.00395 


(0.015, 0.060) 


[30] 


fas 


0.363 ±0.119 


(0.000, 0.85) 


[30] 



Table 2. Nuisance parameters adopted in the scan of the cMSSM parameter space, indicating the mean and 
standard deviation adopted for the Gaussian prior on each of them, as well as the range covered in the scan. 
For scans in which the hadronic and/or astrophysical nuisance parameters have been held fixed, their value 
corresponds to the mean indicated here. 

details): 

Kep)-«, (2.D 

where D is the set of experimental data used to constraint the parameter space, while p(0|D) 
on the l.h.s. of Eq. (2.1) is the posterior probability density function (pdf), p(D\&) is 
the likelihood function (also indicated as £(©)) and p(0) is the prior pdf, representing 
our knowledge of the parameters before the data is taken into account. Bayes' theorem 
expresses how this prior information is updated by means of the experimental data through 
the likelihood function, so that the final pdf for depends both on the priors and on the 
likelihood. The evidence (or model likelihood) p(D) only depends on the data so, for our 
purposes it represents a normalization factor and can thus be dropped. Different choices can 
be made for the prior probability in Eq. (2.1). However if the pdf at the end still exhibits 
a residual dependence on the priors, this represents a clear sign that the experimental data 
employed are not constraining enough to overcome the priors and, as a consequence, the 
resulting pdf should be considered with care [21, 33]. Following Refs. [21] and [22], in order 



to assess the robustness of our inference we consider two sets of priors: one that is uniform in 
the cMSSM masses ("flat" priors) and one that is uniform in their log ("log" priors). Both 
sets of priors are uniform on Aq and tan/3, and the ranges are indicated in Table 1. For the 
SM, astrophysical and hadronic nuisance parameters our priors represents current available 
measurements of those quantities, and are indicated in Table 2 (further details are given in 
Section 3.2 below). Furthermore, we also present results for the profile likelihood, which is 
in principle independent of the prior [22] . 

When presenting results for 1-dimensional or 2-dimensional subsets of the parameter 
space, the remaining parameters can be eliminated in two different ways: 

• By marginalizing over them (i.e., integrating over - Bayesian), resulting in the 1- 
dimensional (or 2-dimensional) posterior pdf for the i-th parameter: 

p(0i|D) = f p(&\r>)d0 1 ...de i „ l dQ i+1 dO n , (2.2) 

• By maximizing over them (i.e., profiling over the likelihood - Frequentist), resulting in 
the 1-dimensional (or 2-dimensional) profile likelihood function: 

£(©,) = max £(&). (2.3) 

01,...,Oi_l,Oi+l]...,OB 

Inferences resulting from those two approaches have, in general, a different meaning and 
may produce different results: the definition of the profile likelihood is more suited to de- 
termine possible small regions in the parameter space with large likelihood values, while, 
on the contrary, the pdf, integrating over hidden directions, correctly accounts for volume 
effects. Usually the largest amount of information on the structure of the scanned space 
derives from the comparison of the two [22, 33]. For both choices of priors, we use a mod- 
ified version of the public SuperBayeS vl.5.1 package [34] 1 to obtain posterior samples of 
Eq. (2.1), adopting MultiNest v2.8 [35, 36] as a scanning algorithm. We use running pa- 
rameters as recommended in Ref. [22], namely a number of live points n\[ ve = 20,000 and 
a tolerance parameter tol = 10 -4 , which are tuned to obtain an accurate evaluation of the 
profile likelihood. Our final inferences for each of the log and flat priors are obtained from 
chains generated with approximately 13 (41) million likelihood evaluations for scans with 
fixed (varying) astrophysical and hadronic nuisance parameters, with an overall efficiency of 
about 5%. The profile likelihood plots are derived from merged chains from both of our log 
and flat priors scans, as advocated in Ref. [22], and are obtained from over 1 million samples. 

2.2 The likelihood function 

The experimental data D included in the likelihood are summarized in Table 3. Observ- 
able quantities already included in Ref. [21] have been used here with the same likelihood 
function as in Ref. [21] (for the ALL setup), with updates to central values and standard 
deviations where applicable as indicated in Table 3. In particular, we include LEP bounds 
on the Higgs mass, precision tests of the SM and an updated constraint on the DM relic den- 
sity from WMAP 7- years data (assuming that neutralinos are the sole constituent of DM). 



For this paper, the public SuperBayeS vl.5.1 code has been modified to interface with with DarkSUSY 
5.0.5 [37, 38], SoftSUSY 2.0.18 [39, 40], MicrOMEGAs 2.0 [41, 42], Superlso 2.4 [43, 44] and SusyBSG 1.4 
[45, 46]. 



We have also made the following modifications with respect to the treatment in Ref. [21]: 
the BR(B — > X s j) is now computed with the numerical code SusyBSG [45], which uses the 
full NLO QCD contributions, including the two-loop calculation of the gluino contributions 
presented in Ref. [47] and the results of Ref. [48] for the remaining non-QCD tan /3-enhanced 
contributions 2 . The other I?(.D)-physics observables summarized in Table 3 have been com- 
puted with the code Superlso (for details on the computation of the observables see Ref. [44] 
and references therein). The anomalous magnetic moment of the muon, <5a^ USY , is also com- 
puted with Superlso. Both Superlso and SusyBSG have been integrated in a modified 
version of SuperBayeS vl . 5 . 1. 

In addition, we also include the recent constraints from SUSY searches at the LHC. 
We follow the analysis presented by the ATLAS collaboration in Ref. [49] , since it is the one 
associated to the more stringent constraints: a data sample corresponding to a luminosity 
of 35 pb _1 for a center-of-mass energy of yfs = 7 GeV is used to search for a SUSY signal 
in events with an isolated electron or muon with high transverse momentum, at least three 
hadronic jets and missing transverse momentum. No significant deviations are found with 
respect to SM predictions, and this is used to derive the constraints on the (?n<o> m i/2) plane 
indicated in Fig. 2 of Ref. [49] as a solid red line. We note that while this 95% exclusion 
line has been obtained assuming values of tan/3 = 3.0 and Aq = GeV, the final result is 
not expected to change much for different values of these two parameters since it is obtained 
from decay channels that are fairly independent on the value of tan/3 and A. We include 
those exclusion limits by assigning a likelihood equal to zero to samples falling below the 
exclusion line. 

The likelihood function for the direct detection constraints is a central ingredient of this 
work, and is presented in more detail in the next section. 

3 Constraints from direct detection data 

3.1 Calculation of the recoil rate 

Before discussing our implementation of direct detection results in the likelihood, we briefly 
review the calculation of the recoil rate, with particular emphasis on the astrophysical un- 
certainties (see also the discussion in Refs. [29, 60, 61]). 

The rate of recoil events produced by DM interactions with a nucleus of mass m^ is 
given by 

§ = ^-[ <W(v + ^M). (3.1) 

ah m x rriN Jy Q-E 

Here, p\ oc is the local DM density, v the DM velocity in the detector rest frame, VEarth the 
Earth velocity in the rest frame of the Galaxy, /(w) the DM velocity distribution function 
and dcr x -N/dE is the differential cross section for the interaction between the neutralino 
and the nucleus. The integral is over the volume V of phase space for which v is smaller 
than the escape velocity v esc and larger than the minimal velocity v m [ n (E) able to produce 



a recoil with energy E: v m - m (E) = Jm^E/2p? N , where [in = m x m^{m x + m^) l is the 
neutralino-nucleus reduced mass. 



2 For certain choices of the input parameters it may happen that the masses of two particles are accidentally 
close to each other. In some of these cases, this leads to numerical instabilities which usually are circumvented 
by compiling the program to quadruple precision. However, we have observed that in rare occasions numerical 
instabilities may persist. 



Observable 


Mean value Uncertainties 


Ref. 




\x a (exper.) 


r (theor.) 




M w [GeV] 


80.398 0.025 


0.015 


[50] 


sin 2 9 e ff 


0.23153 0.00016 


0.00015 


[501 


^USY x 10 io 


29.6 8.1 


2.0 


[511 


BR{B -»• X sl ) x 10 4 


3.55 0.26 


0.30 


[52] 


AM Bs [ps- 1 ] 


17.77 0.12 


2.40 


[531 


BR(B u ^tu) 
BR(B u ->tv)sm 


1.28 0.38 


- 


[52] 


A _ x 10 2 


3.6 2.65 


- 


|54| 


BR(B^Dru) ln2 
BR(B->Dev) x iu 


41.6 12.8 


3.5 


[55] 


-R«23 


1.004 0.007 


- 


[56] 


5i?(L> s -> ri/) x 10 2 


5.38 0.32 


0.2 


[52] 


Si2(D B -> jti/) x 10 3 


5.81 0.43 


0.2 


[52] 


BR(D -> //i/) x 10 4 


3.82 0.33 


0.2 


[52] 


fi x /i 2 


0.1123 0.0035 


10% 


[57] 




Limit (95% CL) 


r (theor.) 


Ref. 


BR(B s ^n+fi-) 


< 5.8 x 10" 8 


14% 


[58] 


m h 


> 114.4 GeV (SM-like Higgs) 


3 GeV 


[59] 


t 2 


f(m h ) (see Ref. [20]) 


negligible 


[59] 


nig 


> 375 GeV 


5% 


[271 


ITlg 


> 289 GeV 


5% 


[271 


other sparticle masses 


As in table 4 of Ref. [20]. 






m ,m 1/2 


LHC exclusion limits, see text 




[49] 


SI 
% ~~ a x-p 


XENON100 exclusion limits, see text 


[24] 



Table 3. Summary of the observables used for the computation of the likelihood function Upper part: 
Observables for which a positive measurement has been made. <5a?° = a™ p — o^ M denotes the discrepancy 
between the experimental value and the SM prediction of the anomalous magnetic moment of the muon (g— 2)^. 
For each quantity we use a likelihood function with mean /j, and standard deviation s = \/a 2 + t 2 , where a is 
the experimental uncertainty and r represents our estimate of the theoretical uncertainty. Lower part: Ob- 
servables for which only limits currently exist. The explicit form of the likelihood function is given in Ref. [20], 
including in particular a smearing out of experimental errors and limits to include an appropriate theoretical 



uncertainty in the observables. nth stands for the light Higgs mass while ( h 
where g stands for the Higgs coupling to the Z and W gauge boson pairs. 



g {hZZ)mssm/g (hZZ) S M, 



We will focus here only on the case of spin independent (SI) interaction, assuming also 
that protons and neutrons are characterized by the same cross section. Therefore: 



do\ 



-N 



dE 



(v,E) 



m N SI 



2/x 2 « 2 



4L P X>^ 2 (^)> 



(3.2) 



where fi p is the neutralino-proton reduced mass. The sum is over all the isotopes present 
in the detector, each one with its own abundance factor Wj. Ai is the atomic mass number 
of the i-th nuclear species and F(Ai,E) the nuclear form factor, which we have taken from 
Ref. [62]. 



3.2 Astrophysical and hadronic nuisance parameters 

Regarding the velocity distribution function entering Eq. (3.1), we account for uncertainties 
in the halo model by considering the following parametrization: 



/(v) = /o 



,2 "i k 



«pI*^I-i 



(3.3) 



The distribution function depends on three parameters, (k,v esc ,Vd), and fo is simply a nor- 
malization factor. Eq. (3.3) has been introduced in Ref. [63] as an Ansatz able to reproduce 
the phase space structure of iV-body simulated DM halos and, compared to the standard 
Maxwellian distribution, it predicts a smaller number of DM particle in the high velocity 
tail, which entails less recoil events. For simplicity we reduce the number of parameters in 
the velocity distribution function to just two, fixing throughout k = 1. 

From Eqs. (3.1) and (3.3), astrophysics affects the number of predicted recoil events 
through the dependence on the three velocities (v Eaith , v csc , Vd) and the local density p\ oc . As 
done in Ref. [29], we note that u Ear th can be reasonably approximated by the value of the local 
circular velocity, v\ ST , so we assume t> Earth = ^i sr . The value of these "astrophysical nuisance 
parameters" is constrained using a Gaussian prior with mean and standard deviation as 
motivated in Ref. [29], and given in the central part of Table 2. In order to derive a prior 
for the velocity dispersion, Vd, we use the relationship Vd = a/3/2-uo and, as in Ref. [29], the 
mean velocity vq is assumed to be equal to v\ ST . However, in our scans we allow both Vd and 
v\ sr to vary as free, independent parameters (in order to be conservative), within the range 
given in Table 2. 

Finally, cr^_ p is proportional to the square of the effective coupling of the neutralino to 
the proton, f p : 

A = E frA + ^fTQ E — • ( 3 - 4 ) 

m„ ^— ' m Q 27 t—' mq 

q=u,d,s q=u,d,s 

where the quantities fx q represent the contributions of the light quarks to the mass of the pro- 
ton, and are defined as fx q = fnp{p\m q qq\p) . The second term corresponds to the interaction 
of the neutralino with the gluon scalar density in the proton, with /tq = 1 — ^2„ =u d s hq- 
The coefficients a s q correspond to the neutralino-quark scalar couplings in the low-energy 
effective Lagrangian and are calculated for each point in our scan. The hadronic matrix 
elements fx q are determined experimentally in nuclear measurements and the uncertainty in 
their value (especially that on frs) has a direct impact on the computed cross section. We 
therefore consider them as nuisance parameters in the scan and we constrain them with an 
informative Gaussian prior, with mean and standard deviation as indicated in Table 2, taken 
from [30] (for an analysis of the capabilities of future ton-scale experiments to constrain dark 
matter including hadronic uncertainties, see [31]). 

3.3 Implementation of XENON100 data 

We have now assembled all the necessary elements to implement the constraints from direct 
detection experiments to our likelihood. The most stringent results to date have been recently 
obtained by the XENON100 experiment, which is located at the Laboratori Nazionali del 
Gran Sasso in Italy. DM is searched for with a target of dual-phase (liquid/gas) xenon in an 
environment of extremely low background. A particle interacting with the detector produces 
a primary scintillation signal SI in the liquid xenon. Subsequently, ionization electrons drift 



towards the region with gaseous xenon producing the secondary scintillation signal 52. The 
background is largely reduced by considering only single scattering events, while the S2/S1 
ratio is employed to discriminate nuclear recoils (produced by backgrounds or candidate 
WIMPs event) from electronic recoils (from 7 or /3 background). 

The latest results [24] have been obtained with an effective volume of 48 kg and 100.9 
days of exposure. The signal region is defined by a series of quality cuts determined blindly, 
in an energy window between 8.4 and 44.6 keV, from the requirement of having from 4 to 
30 photoelectron events (PE). In this region, 3 events have been detected with energies of 
12.1 keV, 30.2 keV and 34.6 keV. The expected number of events from the background is 
b = 1.8 ± 0.6. In our analysis, we adopt a fixed central value for the background. We have 
checked that marginalizing over the background uncertainty has a negligible impact on our 
results. 

In this work, we do not employ the energy distribution of the observed events, but 
only their total number. Having observed N = 3 events and with an expected background 
of b events, the likelihood function for the parameters is a Poisson distribution (up to a 
normalization constant, Lq) 

\N 
£xeioo(©) oc p(N\X) = L -^-e- x , (3.5) 

where A = s(0) + b and the expected signal s(0) is obtained as described below. The 
normalization constant Lq is chosen in such a way that ln£xeioo = when X = N, i.e., when 
the predicted signal plus the expected background match the observed number of counts. 
The probability distribution of the number of PE events, n, is obtained by convolving the 
differential recoil rate with a Poisson distribution with mean S1(E), accounting for the 
acceptance cuts and other quality cuts of the data via a function C(E)~- 



l ^(( E ) P (n\Sl(E))dE , (3.6) 



dR f°° dR 
dn 



where C(-^) is a function parameterizing the acceptance of all data quality cuts, P(n\Sl(E)) 
is a Poisson distribution for n with mean S1(E), and 



>-> 11 



S1(E) = L y ^EL cS (E) (3.7) 



gives the number of PE produced by an event with recoil energy E (we neglect the uncertainty 
due to the finite single-electron resolution of the photomultipliers, which is small, opmt — 0.5 
PE). In the above equation, S m = 0.95 and S ee = 0.58 are the scintillation quenching factors 
for nuclear and electronic recoils, respectively; L y = 2.20 ± 0.09 PE/keV ee is the light yield 
(whose uncertainty we neglect, as it is small) and L e ^{E) is the scintillation efficiency, which 
we take to be the best fit line in Fig. 1 of [24]. We neglect the uncertainty in L e g(E), which 
however is only important for very light WIMP mass, m x < 10 GeV, as argued recently in 
Ref. [64]. As such small values of the neutralino mass are not realized in the cMSSM, our 
assumption has a negligible impact on our results. The expected total number of events is 
obtained by summing over all possible number of PE reaching the photomultipliers, within 
the acceptance window PE m j n = 4 and PE max = 30: 

- £ <£■ < 38) 

n=PE min 





Fixed nuisance 


Incl. nuisance 




Flat priors 


Log priors 


Flat priors 


Log priors 


Input parameters 


m [GeV] 


108.52 


88.68 


220.51 


98.26 


m 1 / 2 [GeV] 


390.96 


370.54 


582.41 


363.28 


A) [GeV] 


-326.06 


-286.55 


990.39 


-278.95 


tan/3 


13.94 


11.10 


24.29 


13.48 


M t [GeV] 


173.383 


173.129 


173.209 


173.617 


m h (m b ) MS [GeV] 


4.210 


4.236 


4.212 


4.271 


^(Mz)^]- 1 


127.956 


127.968 


127.959 


127.961 


a s (M z ) M5 


0.118 


0.119 


0.117 


0.117 


Observables 


m h [GeV] 


115.1 


114.2 


114.7 


114.5 


m x [GeV] 


157.6 


148.5 


238.8 


145.6 


8al USY x 10 10 


18.54 


17.22 


29.20 


20.89 


5i?(5 -> X s ~i) x 10 4 


3.83 


3.36 


3.16 


3.59 


A - x 10 2 


8.78 


8.42 


8.25 


8.45 


BR(D S -► tv) x 10 2 


4.82 


4.82 


4.82 


4.82 


BR(D -»• pi/) x 10 4 


3.86 


3.86 


3.86 


3.86 


°f- P [Pb] 


3.2 x 10" 9 


3.5 x 10" 9 


4.0 x 10" 9 


3.5 x 10" 9 


4-p [Pb] 


8.0 x 10~ 7 


1.1 x 10~ 6 


7.9 x 10~ 7 


1.1 x 10~ 6 


ft x /i 2 


0.1142 


0.1142 


0.1096 


0.1105 



Table 4. Best fit points for the flat and log prior scans. The columns "Fixed nuisance" are for the 
case where astrophysical and hadronic parameters have been fixed to their fiducial values, while "Incl. 
nuisance" denotes the case where they have been included in the scan and profiled over. 



where e is the exposure. We make the further simplification of assuming an energy-independent 
C, and for the combined value of £e we take the effective (post-cuts) exposure of 1481 kg 
days [24]. 

We have checked that the above simplified treatment reproduces fairly accurately the 
exclusion limit obtained in Ref. [24] in the m x — 0"5L„ plane. As we do not make use of 
spectral information, we expect our limit to be more conservative than the one derived using 
the more sophisticated likelihood in Ref. [24, 65]. Indeed, for the same choice of astrophysical 
nuisance parameters as in Ref. [24], we find a 90% limit err} = 0.85 x 10 -8 pb for a WIMP 



- 50 GeV, which is to be compared with the tighter limit found in Ref. [24] , 
"pb. 



mass m 

o-Si_ p = 0.70 x 10 



4 Results 

4.1 Impact of LHC data 

We start by showing in Fig. 1 the impact of the new LHC results on the posterior pdf (for both 
choices of priors) and the profile likelihood of the cMSSM. As one can see by comparing the 
empty grey contours (without LHC data) with the black contours (including LHC limits), 
the impact of current LHC limits is actually fairly modest, since the largest part of the 
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Figure 1. Black contours: posterior pdf (upper panels, for flat and log priors) and profile likelihood (lower 
panels) for the cMSSM parameters, including all present-day constraints (WMAP 7-years and LHC first 
results included), except XENON f00. From the inside out, contours enclose 68%, 95% and 99% of marginal 
posterior probability (top two rows) and the corresponding profiled confidence intervals (bottom panels). The 
black cross represents the best fit, the black dot the posterior mean (for the pdf plots). Parameters describing 
astrophysical and hadronic uncertainties have been fixed to their fiducial values. Blue contours represent the 
constraints obtained without inclusion of LHC data. In the plots on the left, the dashed/green line represents 
the current LHC exclusion limit, while in the right-most plots the red/dashed line is the 90% exclusion limit 
from XENON100, from Ref. [24], rescaled to our fiducial local DM density of p loc = 0.4 GeV/cm 3 . 

favoured region under the remaining data is located at larger values of mo and mi^ (see also 
Refs. [16, 17] for a more detailed analysis). 

The posterior pdf's in this figure exhibit the familiar features of the cMSSM (see e.g. the 
discussion in Ref. [21]). We observe in particular two distinct areas in the mo, irii/2 plane: the 
neutralino-stau coannihilation region, that appears as a narrow band for small scalar masses, 
and a wide area which extends to sizable scalar masses, corresponding to the so-called focus 
point (FP) region [66, 67] (also know as hyperbolic branch [68]). The latter is particularly 
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Fixed nuisance 


Incl. nuisance 


Observable 


Flat priors 


Log priors 


Flat priors 


Log priors 


SM nuisance 


0.204 


0.750 


0.087 


1.247 


Astro nuisance 


N/A 


N/A 


1.060 


0.522 


Hadronic nuisance 


N/A 


N/A 


0.630 


0.667 


M w 


0.986 


1.072 


1.037 


0.828 


sin 2 9 eff 


0.105 


0.900 


0.096 


0.150 


5al USY 


1.755 


2.202 


0.002 


1.089 


BR(B -> X sl ) 


0.527 


0.228 


0.964 


0.008 


AM Ba 


0.274 


0.262 


0.276 


0.315 


BR(B u ^tu) 
BR(B u ^rv) S M 


0.647 


0.614 


0.780 


0.656 


A _ x 10 2 


3.816 


3.306 


3.079 


3.356 


BR{B^Dtu) 
BR(B^Deu) 


0.818 


0.802 


0.843 


0.820 


Rl23 


0.345 


0.340 


0.367 


0.347 


BR(D S -> rv) 


2.224 


2.220 


2.238 


2.225 


BR{D S -> /«/) 


3.109 


3.106 


3.122 


3.109 


£i?(L> -»■ /ii/) 


0.008 


0.084 


0.008 


0.008 


O x /i 2 


0.025 


0.026 


0.056 


0.022 


m/i 


1.022 


1.482 


1.216 


1.297 


XENON100 


0.042 


0.110 


0.0 


0.165 


LHC 


0.0 


0.0 


0.0 


0.0 


Sparticles 


0.0 


0.0 


0.0 


0.0 


Bfl(5 fl -)■ /i+^-) 


0.0 


0.0 


0.0 


0.0 


Total 


15.91 


16.62 


15.86 


16.83 


X 2 /dof 


1.77 


1.85 


1.76 


1.87 



Table 5. Breakdown of the total x 2 by observable for the best fit points in our scans, for both 
cases where the nuisance parameters have been fixed (two left columns) and profiled over (two right 
columns). 



prominent in the case of flat priors (top panels) , due to the "volume effects" associated with 
this prior discussed in Ref. [21] . Along the FP branch the absolute value of the \i parameter is 
of order of the electroweak scale despite scalar particles being very massive (generally several 
TeV), as long as gaugino masses are not too large. The smallness of \x has been interpreted 
as corresponding to low values of the fine-tuning of the model and is due to a "focused" 
behaviour of the Renormalization Group Equation for the Higgs mass parameters. 

This FP region has interesting consequences for the nature of the dark matter. Since \x 
is of the same order as the gaugino mass parameters, the neutralino exhibits a significant Hig- 
gsino fraction (it is a mixed Higgsino-bino state). As a consequence its relic density decreases 
[69], making it possible to satisfy the experimental constraint on the abundance of dark mat- 
ter. The increased Higgsino fraction also leads to a large spin- independent neutralino-nucleon 
cross section, since Higgs exchange along a t-channel becomes more effective (it is propor- 
tional to the product of Higgsino and gaugino components). The cross section can thus be 



as large as a 



si 
x-p 



10 pb for neutralino masses as large as 1 TeV, and in fact the focus 



point region is easily identified in the m x a^_ p plane (see the top right panel of Fig. 1) as a 
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plateau in the region CyL p ~ 10 -7 — 10" 8 pb. 

As we shall see, the latest XENON100 direct detection data set stringent constraints on 
this region. This can already be seen from the right most panels in Fig. 1, which display the 
constraints (without inclusion of direct detection data) in the plane spanned by m x and <7yL p : 
especially for the flat prior case, the FP region appears as a sizable island of high probability 
above the dashed red line, which represents the 90% limit from the most recent XENON00 
data (for standard astrophysical assumptions - we shall relax those in the next section). 
Therefore we expect that the inclusion of the direct detection constraints will further reduce 
the viability of the FP region, as we demonstrate below. 

Turning now to the profile likelihood results (bottom panels in Fig. 1), we see that the 
FP region is excluded at the 95% level, but still viable at the 99% level (outer-most black 
contours), which we are able to map out thanks to our high-resolution scan. Most of the 
favoured region however lies below ~ 1 TeV for both the scalar and gaugino masses. As a 
consequence, we observe in the m x — cr^ (bottom right panel) that the cMSSM parameter 
above the current nominal XENON100 limit is for the most part disfavoured at the 99% level 
from a profile likelihood perspective. We also point out that the 99% region from the profile 
likelihood is much wider than could be assumed just by qualitatively extending either the 
68% or the 95% range, and this owing to the highly non-Gaussian nature of the tails of the 
distribution. Our results therefore indicate that a high-resolution scan is necessary to map 
out the tails of the profile likelihood with sufficient accuracy in order to delimit the 99% 
region, whose extent is much larger than would be inferred by assuming an approximately 
Gaussian distribution from the 68% region. Finally, it is interesting that the extent of the 
99% profile likelihood region is actually very similar to the 99% posterior region from the flat 
prior scan (compare outer-most contours of the middle and bottom panels of Fig. 1). 

We now turn to discuss the best fit points identified in our scans (see Table 4). The 
contribution of each observable to the best fit x 2 is shown in Table 5, where the x 2 is defined 
as 

X 2 ^-2^1n£ 4 , (4.1) 

i 

and L{ is the likelihood for each observable. The likelihood Li is normalized in such a way 
that for Gaussian-distributed data points In Li = when the theoretical value matches the 
observed mean in Table 3. For observables for which only limits exist, the likelihood is 
normalized so that In Li = asymptotically above/below the stated exclusion limit. The one 
exception to this is the XENON100 likelihood, which is normalized as explained above, and 
therefore there is a residual contribution to the x 2 °f ~ 0-57 units when s(&) — > (i.e. for 
very small interaction cross-sections, c^-p <C 10 -9 pb). In evaluating the number of degrees 
of freedom (dof), we only count as "active" Gaussian data points (13 in total, from Table 3), 
for 4 cMSSM free parameters (nuisance parameters are not counted as each one of them is 
independently constrained), thus giving dof = 9. Despite the leakage of probability into the 
FP discussed above, it is interesting that the best fit points for both the flat and log prior 
scans are very well compatible, and lie in the bulk region, in contrast to what had been found 
in Ref. [22, 70], who reported best fits in the FP region based on an older combination of 
data. Despite this, as discussed above the FP region cannot be excluded at the 99% level 
either from a Bayesian or a profile likelihood perspective unless one includes in the likelihood 
direct detection data (see below). 

The observables contributing the most to the total x 2 m Table 5, namely the branching 
ratios of BR(D S — > rv) and BR(D S — > jjlu) are sensitive to new physics mainly through 
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Figure 2. As in Fig. 1, but now black contours include XENON100 data (for fixed hadronic and astro- 
physical uncertainties), while the blue empty contours show for comparison the case where no direct detection 
data are included (from the inside out: 68%, 95% and 99% regions). We observe a strong suppression of the 
viability of the FP region in the Bayesian posterior (top and middle row), and a better agreement between 
the posterior distributions and the profile likelihood. Notice that the XENON 100 90% limit (red/dashed 
line) has been included only to guide the eye, as our implementation of the XENON100 data is slightly more 
conservative than the procedure adopted in Ref. [24]. 

the mass of the charged Higgs bosons H^ and to some extend to tan/3 [71], in particular 
for low m H ± and large tan/3 values. Our best fit points have relatively large m#± and 
low/intermediate tan/3, and therefore the theoretical values for those processes are SM- 
like and have a discrepancy at the 1.5 a level with the experimental values. For the isospin 
asymmetry An- new physics contributions become important to large tan /3 and low mi/2 [72] . 
Hence the predictions for our best fit points are again SM-like, showing a discrepancy at the 
2 a level with the experimental central values. In constrast, relatively good fits can be found 
in the bulk region for both the anomalous magnetic moment of the muon and BR(B — > X s -f), 
as can be seen from their relatively low contribution to the total x 2 - 
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4.2 Impact of XENON100 for fixed nuisance parameters 

We show in Fig. 2 the impact of further adding XENON100 data, and compare the resulting 
constraints with the case where all constraints are included except XENON100 (gray con- 
tours). We keep for the moment all astrophysical and hadronic nuisance parameters fixed to 
their fiducial values, before moving on to the full analysis. The impact of direct detection 
data on the posterior pdf with flat priors is fairly dramatic, as can be seen in the three panels 
of the first row. The most obvious consequence is that the FP becomes strongly disfavoured. 
In the case of log priors (second row of Fig. 2), the impact of direct detection data is less 
dramatic on the posterior, but the parameter space shrinks significantly also in this case, a 
result that clearly demonstrates the potential of direct detection experiments. Overall, we 
observe that the posterior pdf 's for the two choices of priors are now in much better agree- 
ment with each other (despite some residual volume effects, which are evident e.g. in the 
large concentration of probability density for large tan j3 in the middle panel of the first row) , 
as one would expect when the posterior becomes more and more strongly dominated by the 
likelihood and the impact of priors is reduced (see Ref. [73] for another example). The profile 
likelihood 99% region (bottom panels) is also strongly reduced by inclusion of XENON100 
data, and it now extends somewhat in-between the 99% posterior region from the flat and 
the log prior scans. 

Regarding prospects for DM direct detection in the cMSSM, by comparing the right- 
most panels in Fig. 2 we observed that the extension of the favoured region at 99% in the 
m x ,a^_ plane is now almost independent of the statistical approach and the priors used, a 
point we will return to in more detail below when we discuss prospects for cMSSM discovery 
in various observational channels. 

4.3 Influence of hadronic and astrophysical uncertainties. 

We present the effect of including hadronic uncertainties in the top panels of Fig. 3. A 
comparison with the case where the hadronic matrix elements are fixed (grey contours in 
the figure) shows that although a full analysis leads to slightly more conservative results 
by increasing slightly the extent of the favoured regions, the simplified analysis with fixed 
hadronic parameters is accurate enough to draw fairly accurate conclusions on the impact of 
various searches on the cMSSM parameter space. 

Finally, including both sets of nuisance parameters and marginalizing/profiling over 
them is the most robust way of applying direct detection constraints. This case is presented 
in the bottom panels of Fig. 3. Again, we observe a slight broadening of the contours 
for the Bayesian pdf, compared to the case when both sets of nuisance parameters are kept 
fixed, but the overall impact on the favoured region is quite small. Regarding the profile 
likelihood depicted in the bottom-right panel of Fig. 3, the profiled contours for the case 
including nuisance parameters are actually slightly smaller (at 99%) than when the nuisance 
parameters are held fixed. This is a consequence of the slightly better value for the best fit 
X 2 found in our scan including nuisance parameters (cf Table 5), compared with the overall 
best fit when the nuisance parameters are fixed (Ax 2 = 0.05). A detailed analysis reveals 
that the improvement is largely ascribable to the better (g — 2)^ value found in the former 
scan, which is essentially exactly identical to the experimental central value. For the profile 
likelihood, even a very small improvement in the best fit value has an impact on the contours 
as far out in the tails as 99%, as those are defined with respect to the best fit x 2 ■ However, 
given the numerical uncertainties associated with any scan, we can safely conclude that this 
tightening of the contours is a spurious effect and that the extent of the 99% region remains 
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Figure 3. Top row: Impact of marginalizing (profiling) over hadronic uncertainties as nuisance parameters 
in the posterior pdf (profile likelihood, right-most panels). Parameters describing astrophysical uncertainties 
have been fixed to their fiducial values. Bottom row: Impact of further marginalizing/profiling over astrophys- 
ical uncertainties. For comparison, the blue contours are the case where both the hadronic and astrophysics 
nuisance parameters are fixed to their fiducial value. In all panels, all available data have been applied (in- 
cluding LHC and XENON 100). The reach of the future ton-scale XENON IT experiment is also indicated 
(bottom red/dashed line). 



qualitatively the same when including both hadronic and astrophysical nuisance parameters 
in the scan. 



5 Updated prospects for cMSSM discovery 

We now turn to discuss the implications of our results for detection prospects at the LHC, 
and via direct and indirect detection channels. Fig. 4 shows ID posterior distributions and 
profile likelihoods for the lightest Higgs mass, m^, the gluino mass rrig, the neutralino mass 
mo, the spin-independent scattering cross section cr^_ p and the spin-dependent scattering 

cross section of neutralinos off protons &ZZ-P' 

We notice that the ID marginal distributions for the log prior scan (thick solid/red) and 
the ID profile likelihood (dashed/black) are very similar, while the ID distribution from the 
flat prior scan (thin solid/blue) still shows some residual volume effect. This manifests itself 
e.g. in the shift of the bulk of the probability density to larger gluino and neutralino masses. 
However, a robust result of our scans is that the best fit neutralino mass is fairly small (in the 
range ~ 150—250 GeV), the lightest Higgs very light (just above the LEP exclusion limit), and 
the spin-independent scattering cross section a mere factor of ~ 2 below current XENON100 
limits. This therefore puts our best fit point for the cMSSM easily within reach of the next 
generation of direct detection experiments. In particular, the upcoming scaled-up version of 
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Figure 4. ID marginal pdf for flat priors (thin solid/blue), log priors (thick solid/red) and ID profile 
likelihood (dashed/black) for the gluino mass, the lightest Higgs mass, the neutralino mass, the SI and the 
SD scattering cross section (from left to right). Top panels include all data except XENON100, middle panels 
include XENON100 data and bottom panels further marginalize over hadronic and astrophysical uncertainties. 
The best fit point is indicated by the encircled black cross. 

XENON100, XENON1T, is expected to probe by 2015 practically the entire 99% posterior 
and profile likelihood shown above, reaching a sensitivity better than cr?r_ p = 10 -10 pb in a 
mass range extending from 20 to 300 GeV (see e.g. the recent assessment of the XENON1T 
reach in Ref. [74]). 

Interestingly, the FP is a region that presents promising prospects for indirect DM 
searches, especially for what concerns the detection of high-energy neutrinos from DM anni- 
hilations in the Sun (e.g. Ref. [76, 77]), that should be possible with the IceCube neutrino 
telescope if the <7y-» is larger than ~ 3 x 10~ 5 pb (for an observation time of 5 years, and 
including the densely instrumented region DeepCore at the center of IceCube, see e.g. Ref. 
[78] for details). The fact that XENON100 rules out the FP branch has therefore a negative 
impact on the prospects for indirect detection, and further constrains the possibility to probe 
DM in the cMSSM with astrophysical experiments. Upon inspection of the posterior pdf and 
profile likelihood in the right-most panels of Fig. 4, we see that after inclusion of XENON100 
data the region with the highest pdf and profile likelihood values falls below the sensitivity of 
IceCube, although a detection is still statistically possible in the tails of the distributions but 
now more unlikely. Finally, comparison of the middle and bottom rows of Fig. 4 shows that 
inclusion of the astrophysical and hadronic uncertainties does not change the overall conclu- 
sions dramatically, with only mild shifts to the resulting inferences resulting from inclusion 
of the extra nuisance parameters in the scan (bottom row). 
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6 Conclusion 

We have presented in this paper new global fits of the cMSSM, including the most recent 
constraints from the LHC and the XENON100 experiment. Besides the uncertainties on 
Standard Model quantities, our analysis takes into account both astrophysical and hadronic 
uncertainties that enter in the calculation of the rate of WIMP-induced recoils in direct de- 
tection experiment. We have shown that the FP branch of the cMSSM is robustly ruled 
out, even when all uncertainties are taken into account. These results highlight the comple- 
mentarity of direct detection experiments to collider searches. Indeed, direct detection data 
would still be required to rule out the FP even if the LHC exclusion limits are extended to 
a total integrated luminosity of 7 fb _1 at 7 TeV center-of-mass energy [75], which might be 
achieved by the end of 2012. 

Although our results are specific to the cMSSM, we note that the conditions for the 
occurrence of the FP can actually be extended to a general class of unconstrained Supersym- 
metric scenarios. In models where no universality assumptions on the soft mass parameters 
are made, the FP region can be shifted to either larger or lower scalar mass parameters, but 
the general feature of the smallness of the /i parameter is preserved, as are its implications 
for dark matter. In Ref. [76], for example, a low-energy parametrization of the FP was stud- 
ied that could be used to describe this region irrespectively of the high-energy assumptions. 
Although the analysis of these general constructions is beyond the scope of the present work, 
it seems reasonable to extend our conclusions to the general case, given the similarity of the 
predicted direct detection rates. 

While this work was under completion, Refs. [17, 18] performed a similar analysis of the 
impact of LHC and XENON100 data on the parameter space of several theoretical models, 
including the cMSSM. While our results about the viability of the FP regions are qualita- 
tively similar, our analysis differs from that in the above references in that we studied the 
impact of both hadronic and astrophysical uncertainties (both neglected in Ref. [18], while 
Ref. [17] includes only hadronic uncertainties). Furthermore, we adopt a more sophisticated 
statistical framework than Ref. [18] (where the parameter space was explored using random 
scans, rather than the more sophisticated statistical scans employed here) and perform a 
detaild, quantitative comparison between the Bayesian and profile likelihood results (absent 
in Ref. [17]). 
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